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Scientific  Progress 

The  goal  of  this  STIR  project  was  to  initiate  a  research  program  developing  statistical  inference  and  stochastic  simulation 
tools  for  problems  of  diffusion  in  biofluids.  Our  proposed  work  fit  into  a  broad  biophysical  research  effort  on  human  lung 
biology.  This  research  was  conducted  in  the  context  of  a  multidisciplinary  diffusion  research  group  with  members  in  several 
institutions  such  as  Penn  State  University,  UNC-Chapel  Hill,  Tulane  University,  etc.,  and  whose  fields  of  expertise  comprise 
diverse  areas  such  as  Applied  Mathematics,  Probability  and  Statistics,  and  Experimental  Physics.  The  group  has  been 
attracting  funding  from  federal  agencies  such  as  NIH  and  NSF.  This  grant  consisted  of  initial  aspects  of  the  applied  probabilistic 
and  statistical  component  of  the  overall  work.  In  this  context,  one  central  aspect  of  our  research  effort  was  the  use  and  develop 
methods  within  a  biophysical  framework,  as  opposed  to  taking  a  model-free,  phenomenological  approach.  This  means  that  our 
current  and  future  results  are  directly  relevant  for  the  evaluation  of  models  for  microrheology,  such  as  the  GLE,  and  immediately 
interpretable  in  physical  terms.  This  has  been  intended  to  not  only  maximize  the  impact  of  our  work  in  the  broader 
microrheology  research  community,  but  also  to  enhance  interdisciplinary  interaction  within  our  diffusion  group  and  with  respect 
to  data-driven  biophysical  modeling. 


We  focued  on  the  following  action  points.  For  Spectral  Methods,  we  proposed  to 

•  Exploit  traditional  methods  of  semiparametric  spectral  density  estimation  to  estimate  and  establish  the  asymptotic 
subdiffusivity  from  data  for  single  bead  paths. 

•  Devise  spectral  domain  goodness-of-fit  tests  for  relevant  subclasses  of  biophysical  models  for  the  non-asymptotic  dependency 
structures  of  paths,  confirming  or  casting  doubt  on  tra-  ditional  models  such  as  the  GLE. 

For  Heterogeneity,  we  proposed  to 

•  Test  for  heterogeneity  within  and  across  paths  in  the  subdiffusive  regime  by  developing  the  asymptotics  of  parametric  and 
non-parametric  methods. 

•  Test  for  heterogeneity  across  paths  int  he  diffusive  regime  by  using  extensions  to  the  Kalman  framework  and  standard 
spectral  density  estimation. 

•  Test  for  heterogeneity  within  paths  both  within  and  outside  of  the  GLE  framework  by  using  spectral  methods  and  extensions  of 
switching  methods. 

•  Explore  and  identify  alternative  models  which  incorporate  heterogeneity  including  Markov  switching-type  models. 

We  have  accomplished  a  significan  number  of  our  proposed  aims.  In  particular,  we  have 

•  tested  three  competing  detectors  of  subdiffusion 

•  have  been  surprised  by  the  solid  performance  of  the  most  commonly  used,  MSD 

•  have  shown  the  superior  performance  of  wavelet  regression  in  more  real  world  settings. 

This  work  was  done  supporting  an  undergraduate  and  master's  thesis.  In  addition,  we  have  developed  a  computer  package  for 
the  R  language  which  implements  these  methods  along  with  related  simulation  methods.  This  package  is  extendable,  and  we 
are  working  on  incorporation  of  switching  models. 

In  addition,  we  have  mad  significant  progress  on  the  theoretical  development  of  the  MSD  estimator.  This  estimator  was  derived 
intuitively  and  is  the  most  commonly  used  estimator  to  evaluate  subdiffusion  in  microrheology  type  data.  Prof.  Didier  was  able 
to  demonstrate  the  asymptotics  of  such  an  estimator,  which  displays  a  central  limit  type  behavior  in  some  regimes  and 
Rosenblatt  type  limits  in  others.  This  result  demonstrates  that  there  are  both  interesting  applied  problems  in  this  work,  but 
connections  to  interesting  mathematical  results. 
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